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ABSTRACT 

We report the discovery of a new gravitational lens system. This object, 
ULAS J234311.93-005034.0, is the first to be selected by using the new UKIRT 
Infrared Deep Sky Survey (UKIDSS), together with the Sloan Digital Sky Survey 
(SDSS). The ULAS J23431 1.93-005034.0 system contains a quasar at redshift 0.788 
which is doubly imaged, with separation l'/4. The two quasar images have the same 
redshift and similar, though not identical, spectra. The lensing galaxy is detected by 
subtracting point-spread functions from i?-band images taken with the Keck telescope. 
The lensing galaxy can also be detected by subtracting the spectra of the A and B 
images, since more of the galaxy light is likely to be present in the latter. No redshift 
is determined from the galaxy, although the shape of its spectrum suggests a redshift 
of about 0.3. The object's lens status is secure, due to the identification of two objects 
with the same redshift together with a lensing galaxy. Our imaging suggests that the 
lens is found in a cluster environment, in which candidate arc-like structures, that 
require confirmation, are visible in the vicinity. Further discoveries of lenses from the 
UKIDSS survey are likely as part of this programme, due to the depth of UKIDSS 
and its generally good seeing conditions. 
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1 INTRODUCTION 

About 150 gravitational lens systems are now known in 
which a background galaxy or quasar is multiply imaged by 
the action of a foreground galaxy. Gravitational lens systems 
are potentially important because they probe the matter 
distribution of the lensing galaxy independently of whether 
the matter is light emitting (e.g. Cohn et al. 2001; Rusin et 
al. 2002; Rusin & Kochanek 2005; Koopmans et al. 2006; 
Gavazzi et al. 2007). In cases where the background source 
is variable, as in the case of lensed quasars, time delays can 
also potentially be used for determination of the Hubble 
constant according to the method of Refsdal (1964); for re- 
cent reviews see e.g. Courbin (2003), Kochanek & Schechter 
(2004) and Jackson (2007). Lensed quasars can also be used 
for studies of quasar structure via microlensing studies (e.g. 
Agol & Krolik 1999; Ofek & Maoz 2003; Morgan et al. 2007). 

For both statistical purposes, and to locate objects im- 
portant for particular problems, larger samples are needed. 
Currently the largest sample of lenses is the SLACS survey 
(Bolton et al. 2006), which is based on the Sloan Digital Sky 
Survey (SDSS, York et al. 2000) and relies on selecting SDSS 
spectra which contain two different redshift systems at ap- 
proximately the same spatial position. It is mainly sensitive 



to lensed extended galaxies, and has proved to be important 
in tackling the problem of mass distributions in lens galaxies 
(Koopmans et al. 2006). The other large (>20 lens) samples 
include CLASS, a radio-selected survey (Myers et al. 2003, 
Browne et al. 2003) and the SDSS quasar survey (Inada et 
al. 2004, Inada et al. 2005, Oguri et al. 2006). Although indi- 
vidual lenses from these surveys do not always provide good 
constraints on the mass distribution in the lensing galaxy, 
they can be used for Ho determination and also for statis- 
tical studies (e.g. Chae et al. 2003; Ofek, Rix & Maoz 2003; 
Mitchell et al. 2005; Oguri 2007; Oguri et al. 2008). 

As already demonstrated, the Sloan Digital Sky Sur- 
vey (SDSS, York et al. 2000) is potentially a very important 
resource, due to its wide-area sky coverage and consequent 
large number of quasars. The latest catalogue from the SDSS 
(Schneider et al. 2005) contains 77429 quasars. With typical 
lensing rates (e.g. Turner, Ostriker & Gott 1984) this might 
be expected to yield more than 100 gravitational lenses. Un- 
fortunately, many of these lenses will have relatively small 
separations, as high-resolution searches with radio interfer- 
ometers reveal a median separation of between l"and l'.'5, 
and substantial numbers with separations smaller than 1". 
The combination of seeing effects (SDSS imaging has a me- 
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dian PSF width of 1'.'4 FWHM) and weak secondary im- 
ages mean that moderate- to small-separation lenses will 
be difficult to detect directly. Nevertheless, 24 SDSS quasar 
lenses are known (Oguri et al. 2006, Inada et al. 2007) from 
searches at larger separation. At smaller separation attempts 
have been made to overcome the difficulties, including the 
use of colour selection to identify candidates which contain 
light from both quasar and lens galaxy (Ofek et al. 2007) 
and the use of radio catalogues together with optical in- 
formation from SDSS (Jackson & Browne 2007). Here we 
report an alternative approach, using new large-area public 
infrared surveys. 

The UK Infrared Telescope (UK1RT) is currently carry- 
ing out a number of large surveys in multiple (YJHK) infra- 
red bands, some in the areas surveyed by the SDSS, with 
the collective name of the UK1RT Infrared Deep Sky Sur- 
vey (UKIDSS). The UKIDSS project is defined in Lawrence 
et al (2007). UKIDSS uses the UKIRT Wide Field Camera 
(WFCAM; Casali et al, 2007) and a photometric system de- 
scribed in Hewett et al (2006). The pipeline processing and 
science archive are described in Irwin et al (2008, in prepa- 
ration) and Hambly et al (2008). 

Of particular interest here is the Large Area Survey 
(LAS), which will eventually cover 4028 square degrees to 
a depth of K—18.2 in seeing of <l'.'2 and with median see- 
ing below I". This 25-50% improvement over SDSS seeing 
is vital for the discovery of a gravitational lens whose sep- 
aration is just smaller than the SDSS point spread func- 
tion. Moreover, the colour information using the combina- 
tion of optical and infrared bands may also be important 
in prioritising candidates and deciding on their lensing na- 
ture. This paper presents the first results of an attempt to 
use SDSS and UKIDSS together for lens discovery, and fur- 
ther work will extend this to the entire SDSS quasar sam- 
ple. Here we report on the selection of a joint UKIDSS- 
SDSS sample, resulting in discovery of a probable gravita- 
tional lens, ULAS J234311. 93-005034.0. In Section 2 we 
describe the candidate selection; in section 3 we describe 
follow-up observations with Keck imaging and spectroscopy 
of ULAS J234311.93-005034.0, and briefly discuss the re- 
sults in section 4. 



2 SELECTION 

We used the sample of 77429 quasars from the SDSS/Data 
Release 5 (DR5, Schneider et al. 2007); this data release has 
a footprint of approximately 8000 square degrees. Within 5" 
of each object we searched for UKIDSS sources using the H- 
band catalogue from the UKIDSS Large Area Survey Data 
Release 2 (Warren et al. 2007) which has a footprint of 685 
square degrees. 6707 objects are found which have cross- 
matches in the UKIDSS survey, which represent the vast 
majority of SDSS quasars within the UKIDSS DR2 area. 
The UKIDSS DR2 area covers approximately 1 deg either 
side of the equator from right ascensions of approximately 
h -3 h 30 m and 22 h 30 m -24 h , with smaller, mostly equatorial 
fields between 08 h and 16 h . 

UKIDSS //-band images of all of these cross-matched 
objects were then individually inspected by eye. Objects 
were selected for further investigation if they appeared 
to have extensions or secondary structure within approxi- 



mately 2", as this is the range within which selection using 
SDSS alone is likely to become incomplete. In cases of doubt, 
all available UKIDSS wavelength bands were scrutinised to 
check that secondary structure was visible in more than one 
colour. No distinction was made on relative colour of pri- 
mary and secondary, as it is possible that the light of one 
of the images is dominated by light from the lensing galaxy, 
which is expected to be closer to the weaker image. 

From the initial sample of 6707 quasars, 56 candidates 
were identified as possible gravitational lenses. Of these, the 
one known lens (J122608. 0-000602, Inada et al. 2008), is 
recovered as a candidat«0- Two other objects were in the 
list of candidates selected on the basis of SDSS colours by 
Oguri et al. (2007), one of which (J001342.4-002413) has 
already been ruled out by followup observations. 



3 ULAS J234311. 9 005034 

One of the promising candidates identified in the selection 
by eye was ULAS J234311. 93-005034.0, a double source 
with approximate separation of f'.'4. This object has a K- 
band magnitude of 16.5 and an optical "best-PSF" r-band 
magnitude of 20.3, from the UKIDSS and SDSS surveys re- 
spectively. Selected UKIDSS and SDSS (g-band) images are 
shown together in Figure I. The object lies in the FIRST 
20-cm radio survey area (Becker et al. 1995), but no radio 
counterpart is detected at a level of >1 mjy (5a). 

3.1 Keck imaging 

Follow-up imaging and spectroscopy was performed using 
the Keck-I telescope and the Low-Resolution Imaging Spec- 
trograph - Atmospheric Dispersion Compensator (LRIS- 
ADC; Oke et al. 1995) on the night of 2008 January 04. 
Imaging exposures totalling 1140 s were obtained in the R 
filter, and 1500 s in the g filter. The seeing at the time, 
measured from nearby stars, was about 1". Fig. 2 shows the 
Keck images in both filters. In both figures, the object is 
clearly visible as a double source, confirming the impression 
from the UKIDSS survey images. 

The 7i-band image was fitted using the GALFIT soft- 
ware of Peng et al. (2002) using two point spread functions, 
derived from a nearby isolated and unsaturated stellar im- 
age, a Sersic profile, and a flat background to represent the 
sky. The positions and fluxes of the two PSF images were 
allowed to vary, as were the position, flux, position angle, 
effective radius and Sersic index of the galaxy. The resulting 
fits are shown in Fig. 3, and give an image separation of l'/4 
between two point sources differing by about 20% in flux. 
It can be seen, by comparing the raw image and the fitted 
points, that additional flux is required in excess of the point 
sources in order to fit the sources properly. According to the 
fit, this can be provided by a galaxy slightly brighter than 
either image, whose major axis lies close to the line joining 
the two images. Table 1 shows the fitted parameters of the 
galaxy and lensed images. 

1 The footprint of the subsequent release, DR3, of UKIDSS con- 
tains three known lenses, all of which were recovered as candi- 
dates. 
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Figure 1. Top: SDSS g image of the field of 
ULAS J234311. 93-005034.0. The image is 101"on a side 
and shows the object as slightly extended in a direction North 
of East. Bottom: UKIDSS images in four colours (Y, Ji, H and 
K). Each image is 15"on a side. 



Figure 2. Keck images in the g-band (left) and i?-band (right). 
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Figure 3. Zoom-in on the i?-band image of the lensed 
quasars ULAS J234311. 93-005034.0. (Top left) the image of the 
ULAS J234311. 93-005034.0 region; (top right) the fitted PSFs 
for images A and B; (bottom left) the fitted galaxy; and (bot- 
tom right) the residual (image minus the fitted PSF and galaxy 
components). 



A similar attempt was made to fit the g-band image, 
but the fit fails to represent both the lensed images and 
the galaxy. This is likely to be because the galaxy is rela- 
tively weaker, and the fit collapses to a state in which both 
the galaxy and one of the images are at the same position. 
However, if we constrain the position of the centre of the 
galaxy in the g-band to be the same as that for the i?-band, 
we obtain a good fit in which the fitted position angle of 
the galaxy agrees well with that in the i?-band image. Ac- 
cordingly, Table 1 presents this fit to the gr-band image. The 
discrepancy between the two colour bands is likely to be due 
to the large errors resulting from the attempt to fit compo- 
nents separated from each other by distances slightly larger 
than the seeing FWHM. 

Although the quasar is included in the SDSS DR5 
quasar catalogue (Schneider et al. 2007), its SDSS magni- 
tude of i = 20.1, after extinction correction, is fainter than 
the magnitude limit of the SDSS quasar lens survey (SQLS, 
Inada et al. 2003, 2005, 2008; Oguri et al. 2006) and is there- 
fore not part of the statistical SQLS sample. 



3.2 Keck spectroscopy 

Spectra were obtained, also using the Keck-I telescope, on 
the night of 2008 January 11. For the spectroscopy the 
400/3400 grism on the blue arm of LRIS (central wave- 
length 440 nm), together with the 400/8500 grating on the 
red arm (central wavelength 770 nm) and a dichroic cut- 
ting at about 560 nm, were used. The total exposure time 
was 1400 s, and a Of.' 7 slit was used throughout. Data re- 
duction followed standard procedures including bias sub- 
traction with 10 coadded bias frames and flatfielding from 
flat-field frames made using a halogen lamp. The extracted, 
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ij-band 


g-band 


Image A 


Magnitude 


20.60 


20.35 




Position /" 


(0.0,0.0) 


(0.0,0.0) 


Image B 


Magnitude 


20.78 


20.91 




Position /" 


(1.39,0.59) 


(1.21,0.51) 


Galaxy 


Magnitude 


19.94 


21.83 




Position /" 


(0.40,0.20) 


(0.40,0.20)* 




PA 


72 


68 




Axis ratio b/a 


0.43 


0.70 




r e /" 


0'.'94 


1"54 




Sersic index 


1.56 


0.48 



Table 1. Fit parameters for the Keck imaging of 
ULAS J234311. 93-005034.0. All parameters are given with 
respect to the brighter image A, including position offsets in 
arcseconds (positive values indicate increasing RA and Dec). 
Position errors from the fit are nominally 2-3 mas, but are likely 
to be larger, particularly for the galaxy position. The galaxy 
position in the (/-band fit, marked with an asterisk, is fixed with 
respect to image A. Magnitudes are determined in GALFIT 
by the fitting of a PSF to the images, and a Sersic law for the 
galaxy. The photometric scale is calibrated by use of Sloan g,r,i 
values for this object and the transformation equations of Jester 
ct al. (2005). 
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Figure 4. Spectra of images A and B of the 
ULAS 234311.93-005034.0 system. Note that the spectrum of im- 
age A has been displaced upwards by 2xl0~ 16 crgcm~ 2 s~ 1 A~ 1 
for clarity. Cosmic rays and telluric features have been removed. 



sky-subtracted spectra are shown in Fig. 4. The two QSO 
images clearly have the same redshift, which is determined 
to be 0.788±0.001 using the narrow [Om] line at 500.7 nm. 
The spectra are, however, not identical. If we use the fit to 
the i?-band image described above to give the ratio of the 
fluxes of the quasar images, we can recover a spectrum of 
the galaxy alone; this is also shown in Fig. 4. This is obvi- 
ously redder than the quasar images, but there is no obvious 
omission or absorption line. The general shape of the spec- 
trum suggests an upturn due to a 4000A break at around 
5500A, implying a redshift of about 0.3; in this case, deeper 
exposures would be needed to detect the Ca H/K absorption 
lines. 



Line/rest wavelength Image A Image B Ratio (A/B) 



Cm] 190.9 nm (narrow) 7.2 3.3 2.19 

Cm] 190.9 nm (broad) 26.4 20.0 1.32 

Mgn 279.8 nm (narrow) 9.7 5.2 1.86 

Mgn 279.8 nm (broad) 9.5 6.3 1.51 

H/3 486.1 nm 6.4 9.0 0.72 

[Om] 500.7 nm 4.2 6.2 0.68 



Table 2. Line fluxes from images A and B separately, and the 
flux ratio A/B, measured by subtracting a polynomial fit to the 
continuum by eye and fitting Gaussians. All fluxes are in units 
of 10 — 15 ergcm~ 2 s — 1 and are likely to be subject to errors of 
between 5 and 10%. The broad component of Mgn in particular 
probably contains large amounts of contamination from Fell. 



4 DISCUSSION AND CONCLUSIONS 

The detection of two images with the same redshift, to- 
gether with clear evidence for a lensing galaxy, convinces 
us that ULAS J234311. 93-005034.0 is a lens system. Fur- 
ther circumstantial evidence can be obtained by considering 
the expected brightness of the lens galaxy from the image 
separation and the Faber- Jackson relation (Faber &: Jackson 
1976). If we use the lens separation together with the Faber- 
Jackson relation as calibrated for a sample of gravitational 
lenses by Rusin et al. (2003), and assume a lens redshift of 
0.32, we obtain an expected .R-band magnitude of 19.7, very 
close to the actual value. 

Nevertheless, it is obvious that the spectra (Fig. 4) are 
not identical, despite the very similar redshift. The contin- 
uum shape is different, but as we have seen it is plausible 
that this is due to different levels of contamination from the 
lensing galaxy. A complicating factor, however, is that there 
may be differential reddening between images A and B. This 
can be deduced from ratios of the broad emission line fluxes 
between image A and B, which decrease systematically with 
wavelength (Table 2). If this reddening takes place in the 
lensing galaxy, and if its redshift is about 0.3, this implies 
an associated E(B— V) of about 0.7 for a Galactic extinction 
law. 

Assuming that the reddening occurs in the lensing 
galaxy, it is then puzzling why the fit to the Keck images 
places the lensing galaxy closer to A than to B. There are 
precedents for this (e.g. CLASS B0218+357; Jackson, Xan- 
thopoulos & Browne 2000, HE1104-185; Wisotzki et al. 
1993) but an alternative explanation is that the fitted cen- 
troid of the galaxy is in error. HST imaging would be useful 
to resolve this matter fully. Reddening of the quasar con- 
tinuum is difficult to assess, due to contamination of both 
A and B spectra by light from the lensing galaxy. An addi- 
tional complication is the possibility of microlensing in the 
lensing galaxy, which may affect both the quasar continua 
and, if the broad line region is small enough, the broad line 
fluxes. It is becoming increasingly clear (e.g. Kochanek et 
al. 2007, Morgan et al. 2006, Keeton et al. 2006, Gaynullina 
et al. 2005) that microlensing is important in understanding 
relative fluxes in images of lensed quasars. 

It is possible to make a lens model to describe the 
system, although it suffers from a lack of constraints. The 
available constraints are the image positions and fluxes (as- 
suming we ignore effects of reddening or microlensing), and 
the galaxy position is known. We assume that positions 



The gravitational lens ULAS J234311. 93- 005034.0 5 




-1 o 1 



Figure 5. Straw-man lens model for the 
ULAS J234311. 93-005034.0 system. The scale on x and y 
axes is in arcseconds. Images A and B are marked, together with 
the centre of the galaxy (G) and the model source position (S). 
The greyscale shows the predicted images and the contours are 
time-delay contours. The measured positions of the images (at 
either end) and lensing galaxy (second from right) are shown, 
as is the inferred position of the source. The model includes 
external shear of magnitude 0.34 and direction —18°, and has a 
fixed ellipticity of 0.5 in PA 72°. The predicted Einstein radius 
of the lens galaxy is 0"58. 

are known to 10 mas and fluxes to 10%, and use a singu- 
lar isothermal ellipsoid to represent the mass profile of the 
galaxy. Fixing the ellipticity at 0.5 and the position angle at 
the same as that of the light distribution, but varying the 
lens Einstein radius, gives x 2 /d.o.f .=595/2. If in addition 
we allow the ellipticity to vary, we obtain x 2 /d.o./.=48/l at 
the cost of an unlikely solution of an almost circular galaxy. 
If we once more fix the ellipticity at 0.5, but allow two ex- 
tra free parameters to represent external shear, we remove 
all the degrees of freedom and therefore obtain an exact fit 
(Fig. 5). This fit has a very large shear magnitude (0.34) 
at a position angle of —18°, and predicts a time delay of 
18 days for Hq=70 kms _1 Mpc _1 , with the A image varying 
first. Interestingly, the predicted shear direction is within a 
few degrees of the direction of the nearest object, a galaxy 
about 10" away and slightly west of south. 

That there may be an effect of the environment on the 
lens model is consistent with the appearance of the Keck 
images. The ULAS J234311. 93-005034.0 system appears to 
lie at the northern end of a cluster of galaxies, the brightest 
one of which lies about 49" South. The redshift of this bright 
galaxy appears in the SDSS database as 0.316, so it is at least 
plausible that the ULAS J234311. 93-005034.0 lens galaxy 
redshift lies at around this value. There are a a number of 
arc-like structures in this cluster, suggestive of weak lensing 
of background objects (see Fig. 2). The velocity dispersion 
of the cluster implied by the presence of a shear of 0.34 at 
the position of the gravitational lens is large (~1400 kms -1 
assuming an isothermal profile for the cluster) but detailed 
modelling with additional constraints is needed in order to 
make a serious estimate of the cluster mass. 

The discovery of ULAS J234311. 93-005034.0 suggests 



that the use of large-area surveys with good spatial reso- 
lution is a good way to discover lenses in existing surveys 
of quasars such as the SDSS quasar sample. The increasing 
coverage of the UKIDSS survey means that ultimately be- 
tween 60000-70000 high-resolution images of quasars should 
be available. For typical angular separation distributions of 
existing lens samples, this should allow the existing sample 
of 24 SDSS quasar lenses to be expanded by at least a factor 
2, and possibly more. 
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